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The Chicxulub impact is thought to have triggered a global winter at the 


Cretaceous-Palaeogene (K-Pg) boundary 66 million years ago. Yet the 
climatic consequences of the various debris injected into the atmosphere 
following the Chicxulub impact remain unclear, and the exact killing 
mechanisms of the K-Pg mass extinction remain poorly constrained. Here we 
present palaeoclimate simulations based on sedimentological constraints 
from an expanded terrestrial K-Pg boundary deposit in North Dakota, United 
States, to evaluate the relative and combined effects of impact-generated 
silicate dust and sulfur, as well as soot from wildfires, on the post-impact 
climate. The measured volumetric size distribution of silicate dust suggests 
a larger contribution of fine dust (-0.8-8.0 um) than previously appreciated. 
Our simulations of the atmospheric injection of sucha plume of 
micrometre-sized silicate dust suggest a long atmospheric lifetime of 15yr, 
contributing toa global-average surface temperature falling by as much as 
15°C. Simulated changes in photosynthetic active solar radiation support 

a dust-induced photosynthetic shut-down for almost 2 yr post-impact. We 
suggest that, together with additional cooling contributions from soot 

and sulfur, this is consistent with the catastrophic collapse of primary 
productivity in the aftermath of the Chicxulub impact. 


The Chicxulub asteroid impact event ~66 million years ago showcases 
aunique opportunity to examine the rate, magnitude and mechanisms 
of extreme and abrupt climate change in Earth’s history’ *. The 45-60° 
inclined impact of a10-to15-km-sized carbonaceous chondrite onthe 
Yucatan Peninsula in Mexico” triggered a chain reaction of events 
ultimately responsible for the Cretaceous-Palaeogene (K-Pg) mass 
extinction and the demise of ~75% of species, including the iconic 
non-avian dinosaurs”. Early studies have suggested that ejecta release 
of crater material into the higher atmosphere”*°, as derived from the 
shock-vaporization of the Yucatan target stratigraphy consisting ofa 


carbonate and evaporite platform on top of a granitic basement” 
(Fig. 1), triggered a global impact winter, leading to a dramatic decrease 
in global annual mean surface temperature by more than 10 °C in the 
first year(s) after the impact?*>? 1°, 

More-recent Chicxulub palaeoclimate modelling studies?" 
investigated with increased detail the short-term climatic effects of 
atmospheric injections of fine-grained ejecta, which we define as a 
group of impact-generated material ejected out of the impact struc- 
ture and characterized by a dominant grain size of generally <10 um. 
For the Chicxulub case, fine-grained ejecta consist largely of three 
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Fig. 1 | Conceptual model of the Chicxulub impact plume showing different 
stages of production, transport and deposition of coarse and fine-grained 
impact-generated ejecta (not to scale). a, Chicxulub impact site, Yucatán 
Peninsula, Mexico. The excavation stage shows a 45° impact of a 10- to 15-km- 
sized bolide from the northeast’. The Yucatán target rocks are displaced to 
forma transient cavity with an uplifted rim and are ejected within an impact 
plume and ejecta curtain”. Glassy impact melt spherules*ć and shocked 
minerals from the granitic basement’ are rapidly expelled out of the crater. b, At 
the terrestrial Tanis K-Pg site“ in North Dakota (United States), -3,000 km north 
of Chicxulub, rapidly moving impact-induced shock waves trigger a seiche- 
inundation event on top of a riverine point bar deposit. Meanwhile, the melt 
spherules are—largely ballistically—emplaced in the event deposit, followed 

by the arrival of shocked mineral clasts, constraining the time of deposition at 
<2 hafter impact”*. c, At the impact site, in the minutes after the excavation 
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stage, the cavity rim collapses whereas rocks in the central area are first uplifted 
and then subside downwards and outwards to form a peak ring®. Meanwhile, 
the impact plume with fine-grained ejecta (silicate dust, soot and sulfate 
particles from S-bearing gases) grows, expands and rises into the stratosphere, 
thereby distributing material around the globe”. d, At Tanis, in the first years 
post-impact, atmospheric settling is taking over with Ni-rich spinel-bearing 
spherules, soot particles derived from wildfires and a large fraction of silicate 
dust“. The panel shows a schematic stratigraphy of a silt-to-sandstone interval 
with glassy spherules, followed by claystone with shocked minerals, Ni-rich 
spinels and fine dust rich in iridium, overlain by a lignite”. The claystone 
interval just below the Palaeogene coal is the main focus in this study*””. 
Panels a,c adapted with permission from ref. 10, Springer Nature Limited. 
Credit: tree icon, Flaticon.com. 


groups (Fig. Ic): (1) sulfur-bearing particles produced by water and 
sulfur-bearing gases, resulting from shock-vaporized evaporites*”®, 
(2) soot particles, probably generated from burning of organic-rich 
target rocks and release from global wildfires?” and (3) silicate dust 
particles, derived from pulverization of the Yucatan crystalline base- 
ment”. Besides these fine-grained ejecta, the collision also released 
other products from the Yucatan target area into the atmosphere, 
such as carbon dioxide, water vapour and methane”. Elevated con- 
centrations of these climate-active gases resulted, after the initial 
impact winter”™, in global warming*’°“; CO, also played a key role in 
post-impact ocean acidification and primary productivity”. However, 
these products were acting predominantly ona much longer timescale 


(in the order of thousands to tens of thousands of years’°”°) than con- 
sidered in this palaeoclimate study. 

The emitted mass and particle properties of impact-generated sul- 
fur, soot and dust remain relatively poorly constrained to date. There- 
fore, many uncertainties remain surrounding the atmospheric lifetime 
of fine-grained ejecta and their short-term climatic effects, which in 
part result from a limited amount of observational data from K-Pg 
boundaries incorporated in palaeoclimate models. Moreover, no com- 
bined palaeoclimate scenario emitting concurrently all fine-grained 
ejecta components has been considered so far for the Chicxulub case”. 
Quantifying relative and combined roles of these fine-grained ejecta 
on the global K-Pg climate crisis is paramount to better understand 
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the severity and duration of the K-Pg impact winter and to identify the 
killing mechanisms responsible for the mass extinction. 

The K-Pg impact winter hypothesis, proposed originally in19807, 
was centred around dust-sized material ejected from the crater into 
the stratosphere and subsequently spread around the globe’. However, 
this dust scenario was rejected in the early 2000s because the fraction 
and mass of this clastic debris at the K-Pg boundary was estimated to 
be too modest to cause an impact winter”. Alternative studies pointed 
out that the post-impact sulfur release constituted a more important 
driver behind the prolonged K-Pg impact winter”, as sulfur-bearing 
gas would have resided in the stratosphere for a period of decades 
in the form of sulfate particles®***”’. Previous studies’*” reported a 
longer cessation of solar irradiation as a result of stratospheric sul- 
fate particles relative to the effect of dust particles, emphasizing the 
strong focus of previous modelling efforts on the role of sulfur release. 
Since the mid-1980s, a soot-driven K-Pg impact winter has also been 
postulated”°. The impact-generated soot particles could also have 
played a dominant role in the global blockage of solar irradiance and 
the prolonged post-impact cooling”, as fine soot constitutes a strong 
sunlight absorber”. Soot and charcoal remains found in the K-Pg 
boundary intervals around the world suggest widespread wildfires in 
the aftermath of the Chicxulub impact”. These wildfires were prob- 
ably formed proximally during the first ‘fireball stage’ and poten- 
tially also elsewhere around the globe due to the thermal radiation 
from the atmospheric re-entry of hypervelocity ejecta®’**. A recent 
molecular-burn marker study revealed that up to 17% of the soot release 
was derived directly from Yucatan target rock, while the remainder 
was sourced from delayed wildfires”. However, the overall extent and 
intensity of these wildfires are still debated”. 

The third type of fine-grained ejecta represents the silicate dust, 
which remains a poorly understood group mainly because of a limited 
number of studies measuring their particle-size distribution. This dust 
is clearly finer than other silicate ejecta such as coarse glassy impact 
spherules (microtektite-like; with diameters ranging from ~300 pm 
to 2mm (ref. 36)) and ejected mineral clasts (for example, shocked 
quartz, zircon, feldspar; diameters ranging from ~30 to 600 pm 
(refs. 9,37)), both of which are associated with ejecta curtain processes” 
(Fig. 1a,b), as well as crystalline impact spherules (microkrystites with, 
for example, Ni-rich spinels; diameters ranging from ~100 to 300 um 
(refs. 36,38)), which are linked to condensation out of the vapour-rich 
impact plume?®?’ (Fig. 1c,d). Previous studies considered two sub- 
groups of silicate K-Pg dust”. The first subgroup consists of material 
with a median diameter of 0.5 um (ref. 39) and is known as submicron 
clastic dust. However, the amount of this subgroup was estimated 
to be <6 x 10" g and was therefore considered less important to the 
post-impact climate change”””’. The other subgroup represents nano- 
particles, potentially iron-rich*®, with a median diameter of 20 nm 
(ref. 39). Their total mass was estimated to be much larger, at 2 x 10° g 
(ref. 39) as an upper limit, and recent climate model results using bal- 
listic ejection of this subgroup led to their rapid atmospheric removal 
within nearly 2 yr after impact”. 

To understand the K-Pg dust conundrum, this study includes 
high-resolution sedimentological constraints from an expanded and 
well-preserved K-Pg sequence from the US western interior (Tanis, 
North Dakota; Fig. 1b,d and Extended Data Fig. 1). In contrast to other 
intermediate or distal K-Pg boundary sites, in which the K-Pg clay 
interval is around 1 cm thick and different types of impact ejecta are 
often mixed, the unique Tanis K-Pg record is expanded and allows for 
subselection of specific time frames during the Chicxulub impact 
ejecta deposition. The ~1.3-m-thick ejecta-bearing sediment pack- 
age emplaced by a Chicxulub impact-generated tsunami-like surge 
(Fig. 1b and Extended Data Fig. 1) was deposited exclusively during 
the period of coarse ejecta accretion, probably corresponding to <2 h 
post-impact*’. Shocked quartz grains have been found only in the 
un-reworked K-Pg claystone that conformably overlies the impact 
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Fig. 2 | Volume-weighted K-Pg boundary grain-size data. The laser-diffraction 
grain-size distribution curve is shown of the silicate lithogenic fraction for the 
uppermost K-Pg boundary claystone level from the Tanis site, just below the 
Palaeocene lignite (see Extended Data Fig. 1 for a detailed stratigraphy). The 
orange curve displays the volume-based grain-size distribution, with the orange 
bar and blue representing -67 and -95% of the data, respectively (corresponding 
to grain-size ranges of -0.8-8.0 um and -0.2-35.0 um, respectively). The black 
dashed line highlights the median grain-size value after a conversion into 

a number density spectrum (see Methods and Extended Data Fig. 2 for the 
converted number density spectrum). 


spherule-bearing Tanis deposit (Fig. 1b,d). On the basis of ballistic 
trajectory simulations, these shocked minerals are expected to arrive 
at Tanis in -2 h after impact’, corroborating the temporal control. A 
positive iridium anomaly exists in this K-Pg claystone (3.8 ppb (ref. 41)), 
whichis linked to impactor material and is comparable in concentration 
to other globally distributed K-Pg localities‘. The uppermost part of the 
K-Pg claystone interval, just below the Palaeogene lignite, is considered 
representative of the final atmospheric settling of the fine silicate dust, 
richiniridium (Fig. 1d). Laser-diffraction grain-size analysis performed 
onthe different strata at the Tanis K-Pg site quantifies the particle-size 
distributions (Methods and Extended Data Fig. 1). This uppermost K-Pg 
claystone interval represents a distinct, uniform and very fine distribu- 
tion witha volume-based median grain size of 2.88 1m (Fig. 2), and after 
converting to anumber density spectrum with a median grain size of 
0.125 um (mass to number density spectrum conversion in Methods; 
Extended Data Fig. 2). These size distribution data are used here as an 
input parameter to simulate the silicate dust scenario. 

Previous atmospheric modelling studies investigating the radia- 
tive effect of silicate dust following the Chicxulub impact event!®?? 
have used either nanometric-sized particles or coarse impact spherule 
data as the bulk fraction of silicate dust injection, both of which are 
estimated to be onthe order of 2 x 10" gas an upper limit. The dominant 
grain sizes between ~0.8 and 8.0 pm measured at the Tanis K-Pg site 
result in substantially slower deposition than those used in previous 
palaeoclimate studies'*”’ (Fig. 3). In this Article, we incorporate this 
sedimentological data in our general circulation model (GCM)”, which 
makes use of the latest Cretaceous palaeogeographic reconstructions” 
(Extended Data Fig. 1a) and the latest insights onthe moment of impact, 
which is suggested to be the boreal spring season on the basis of recent 
data from the Tanis K-Pg site” (Methods). In our palaeoclimate GCM 
simulations (see Methods for details), we injected into the atmos- 
phere the same amount of silicate dust (2 x 10" g (ref. 39)), fine soot 
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Fig. 3 |Number-weighted K-Pg boundary grain-size data in relation to 
atmospheric settling processes. Phase diagram shows grain size versus 

the variation of dry deposition velocity over land (green solid curves) and 

ocean (blue solid curves) for two typical friction velocity values, u., on Earth, 
indicating the intensity of surface wind shear (Methods). The nonlinear 
responses of deposition rates along with particle sizes are controlled by the 
relative contributions of gravitational and aerodynamic forcing, besides the 
Brownian diffusion, impaction and interception of particles (Methods). Previous 
modelling studies!” have used three groups of ejected dust in different sizes. 
The first group is nanometric in size (median of 0.02 um; brown dashed line), 
following a fast Brownian deposition. The second group comprises ejected 
(shocked) quartz grains (defined by ref. 27 as clastic debris; red dashed line) 
witha mean size of 50 um (ref. 37). Together with the third group, the type 
2spherules (microkrystites”) with a median size of 250 um (purple dashed 

line), this follows a fast gravitational settling. These literature data contrast the 
silicate dust fraction measured at the Tanis K-Pg site (number-based median of 
0.125 um with two-thirds of the particles by number density falling in the range of 
~0.08-0.19 um and 95% of the data falling in the range of -0.05-0.30 um; orange 
and blue dashed line and bar). This implies a considerably slower deposition for 
the impact-generated dust, by a factor of -3-4 compared with the nanoparticles 
and bya factor of -10 or at least 100 compared with the ejected quartz grains and 
type 2 spherules, respectively. 


(2.08 x 10" g (ref. 18), using a median diameter of 0.22 um (ref. 39); 
Extended Data Fig. 3) and sulfur (3.25 x 10” g (ref. 16), equivalent to 
6.5 x 10” g SO, (ref. 18)) as used in recent K-Pg modelling studies to 
directly compare the simulation results from our dust model with 
previous works'°*"6°, Using these input parameters, the individual 
scenarios are modelled for silicate dust, soot and sulfur, accompa- 
nied by acombined scenario. In this Article, we carefully document 
the short-term global environmental consequences of the Chicxulub 
impact witha focus on changes in the photosynthetic active radiation 
(PAR; Methods). This solar radiation from 400 to 700 nm is being used 
by organisms inthe process of photosynthesis and therefore represents 
an efficient indication for the state of the biosphere. 


Cooling in the aftermath of the Chicxulub impact 

Our palaeoclimate simulations show global temperatures changing 
seasonally between ~15 and 19 °C before the impact event, in line with 
proxy-based latest Cretaceous temperature reconstructions***, 
despite a larger meridional temperature gradient compared with the 


proxy temperatures (Extended Data Fig. 4). The K-Pg impact winter 
catastrophically disrupted this greenhouse climate, leading to plum- 
meting temperatures in both the marine and especially the terrestrial 
realm (Extended Data Fig. 5). The global-average surface temperature 
dropped by as muchas 25 °C (Fig. 4a). With varying magnitude and 
timescales, each fine-grained ejecta type results in a prolonged global 
cooling inthe initial aftermath of the impact. The sulfur scenario causes 
the lowest surface temperatures, with a global average of -3 °C, fol- 
lowed by the silicate dust scenario, with a global average of ~7 °C. Our 
simulations indicate that both sulfur and silicate dust have a prominent 
influence on the surface thermal forcing, leading to a prolonged impact 
winter up to ~20 yr, in which the first 5-8 yr are the most severe. 


Relative radiative contributions of dust, sulfur 
and soot 

The ejecta type and its atmospheric mass determines the magni- 
tude and timescale of the impact winter (Fig. 4b), related to different 
responses in net short-wave/long-wave radiation (Extended Data Fig. 6 
and Methods) and downward long-wave radiation (Fig. 4c) on the 
surface. The latter indicates the reflected infrared radiation due to 
atmospheric particles and clouds, which may either diminish or further 
enhance surface warming. Indeed, it plays a crucial role in Earth’s cli- 
mate asa vital factor controlling surface warming. The mass of sulfate 
particles in the atmosphere decreases much faster than dust, such that 
the global column-integrated mass of sulfur is reduced to about two 
orders of magnitude lower than that of dust at 3 yr after impact (Fig. 4b). 
Because of the radiatively transparent optical feature of sulfate par- 
ticles”, the sulfur emission is not able to enhance the downward 
long-wave radiation. Injust a few years (Fig. 4c), it drops by 31.4% from 
its pre-impact level, thus allowing more long-wave radiation and caus- 
ing the most cooling of all the impact scenarios (Fig. 4a). The decline 
or enhancement of the downward long-wave radiation is what makes 
the distinction between sulfur and dust cases, respectively, linked to 
the differences in particle optical properties, that is, specific extinc- 
tion, single-scattering albedo and asymmetry factor (Methods). The 
dust-induced downward long-wave radiation shows a rapid increase 
by 64.4% from the pre-impact level to 372 W m?” ina week after impact 
(Fig. 4c), leading to relatively warm surface temperatures more than 
~10 °C compared with the sulfur case (Fig. 4a). Following the initial 3 yr 
after impact, the downward long-wave radiation increases gradually to 
the pre-impact level (210 W m°), lasting below pre-impact level -20 yr 
in all cases except the soot-only case, which reached the pre-impact 
level in ~8 yr (Fig. 4c). This trend is consistent with the gradual global 
warming pattern at the same time frame (Fig. 4a). Our simulation 
results depict that the radiative feedback of the combined emissions 
scenario (Fig. 4c,d and Extended Data Fig. 6), hence the surface cooling 
(Fig. 4a), does not result in the sum of individual ejecta groups but 
rather falls between sulfur and dust cases associated with distinct radia- 
tion responses. It is governed largely by nonlinear interaction processes 
inthe atmosphere, such as the complex dynamics of aerosol transport 
and deposition mechanisms interacting with turbulent atmospheric 
circulations and the water cycle. 


Photosynthetic active radiation cessation 
following impact 

The global-average PAR (Methods) flux inthe sunlight spectrum varies 
seasonally between 80 and 95 W m° in the latest Cretaceous period 
(Fig. 4d). The global mean corresponds to 89 W m”, which is larger 
than the present-day global PAR flux (66-74 W m°; Methods). Because 
of the seasonality, the PAR flux is as high as 160 W m” in the Northern 
Hemisphere just before the impact winter (Extended Data Fig. 7 and 
Methods). The first-day response (Fig. 5) displays the trace of initial 
ejecta expansion through the atmosphere, marking a regional PAR 
deficit onthe surface. Given that the dust and sulfur are ejected from the 
impact target, their traces are initially centralized around the Yucatan 
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Fig. 4 | Temporal evolution of the Chicxulub impact-generated global 
climatic responses. The time evolution is shown from 2 yr of the latest 
Cretaceous towards 25 yr of post-impact conditions, for the individual 
silicate dust, sulfur, soot and combined scenarios. a, Global-average surface 
temperature. b, Global column-integrated fine-grained ejecta mass. 


c, Global-average surface downward long-wave radiation flux. d, Global-average 
PAR flux. On the x axis, the year O refers to the start of the year when the impact 
event occurs. The solid purple dashed line denotes the moment of Chicxulub 
impact, that is, during the boreal spring season”. 


area, unlike the global spread of soot released from wildfires (Methods). 
Within two weeks after impact, the planet undergoes a global shut-down 
in photosynthetic activity in all ejecta scenarios (Extended Data Fig. 8). 
Consistent with global-average responses (Fig. 4d), the massive dust 
ejecta indicates the largest influence on the land-ocean photosynthetic 
activity, leading to an entire PAR suppression for ~1.7 yr (Fig. 5a). Even 
2.2 yr after impact, during boreal summer, the photosynthetic activity 
remains low although in partial recovery with a moderate PAR flux of 
100 W m?” (Fig. 5a), reaching pre-impact values in ~4 yr after impact 


(Fig. 4d). Following the initial suppression of six months (Fig. 4d), the 
PAR deprivation in sulfur and soot cases recovers faster (>1 yr) thanin 
the dust case (Fig. 5a,b,c). This relatively rapid PAR recovery clearly 
documents that, alone, these two factors (sulfur and soot) are much 
less lethal than dust. Moreover, the combined emissions scenario indi- 
cates a global PAR cessation even at ~1.7 yr after impact (Fig. 5d). It 
fully recovers only after ~4 yr after impact, similar to the dust scenario 
(Fig. 4d), implying that the dust is the most lethal for the planet’s photo- 
synthetic cessation. 
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Fig. 5 | Temporal evolution of PAR flux reconstructions after the Chicxulub 
impact. a-d, Land and ocean PAR flux from one day before (pre-impact state) 
and after (post-impact state) in Boreal spring instantaneously towards 5-6 yr 
(annual mean) in the aftermath of the Chicxulub impact for silicate dust (a), 
sulfur (b), soot (c) and combined (d) scenarios, displayed ona latest Cretaceous 


palaeogeographic map (Extended Data Fig. la). The range of the green-white 
colour bar represents the photosynthetically high and low radiative flux. 
Basemaps are based on the latest Cretaceous palaeogeographic data for the 
latest Cretaceous”. 


Implications for the K-Pg climate and biota 

In contrast to previous work*"*”””°, our palaeoclimate simulations, 
relying on robust sedimentological field data at the K-Pg, revealed that 
a massive impact-generated dust plume plays a large part in driving the 
K-Pg climate and biotic crisis. This dust cloud is composed of a large 
number of climate-active particles with grain sizes primarily in the 
range of 0.8-8.0 um by volume (Fig. 2). This information is deciphered 
from laser-diffraction grain-size analysis from the uppermost K-Pg 
boundary claystone of the terrestrial Tanis site (median of 2.88 um), 
which shows a comparable grain-size distribution as the deep marine 
ODP 1212 site in the Pacific Ocean (2-4 um (ref. 46)). Our palaeoclimate 
simulations indicate that this micrometric grain-size pattern results 
in an atmospheric residence time that is much longer (>2 times) com- 
pared with previously estimated nano-sized or coarse particles'*””” 
(Fig. 3 and Extended Data Fig. 9). Among simulated ejecta types, our 
simulations demonstrate a post-impact world where the silicate dust 
represents the largest ejected mass in the initial aftermath of impact 
(Fig. 4b). The sulfur scenario exhibits the fastest deposition with an 


atmospheric lifetime of only -8.5 yr. By contrast, both soot and silicate 
dust resided in the atmosphere considerably longer, for -15 yr after 
impact. Similarly, the combined emissions scenario shows an atmos- 
pheric lifetime of -15 yr, implying a plausible condition for the Chicxu- 
lub impact. The complete recovery from the impact winter took even 
longer, with pre-impact temperature conditions returning only after 
~20 yr (Fig. 4a and Extended Data Fig. 10). This timescale is consistent 
with the recent global iridium layer observations’ from the Chicxulub 
impact structure, where the final atmospheric settling of fine-grained 
impactor material in the dust cloud was estimated to be <20 yr. 

We find that the global darkness and prolonged loss inthe planet’s 
photosynthetic activity occur only in the silicate dust scenario, up to 
nearly 1.7 yr (620 days) after impact (Figs. 4d and 5a-c). This consti- 
tutes a sufficiently long timescale to pose severe challenges for both 
terrestrial and marine habitats. Biotic groups that were not adapted 
tosurvive the dark, cold and food-deprived conditions for almost two 
years would have experienced massive extinctions. This matches the 
palaeontological records, which show that fauna and flora that could 
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enter a dormant phase (for example, through seeds, cysts or hiber- 
nation in burrows) and were able to adapt to a generalistic lifestyle, 
not dependent on one particular food source (for example, deposit 
feeders), generally better survived the K-Pg event”. In addition, our 
combined emissions scenario shows that the photosynthetic recovery 
to the pre-impact levels first occurred in the austral summer season, 
~1.7 yr after impact. This would imply an earlier recovery of the primary 
productivity inthe Southern Hemisphere, consistent with palaeonto- 
logical evidence suggesting lower extinction levels on the Southern 
Hemisphere“, resulting in geographic heterogeneity in extinction 
and recovery from the end-Cretaceous catastrophe. However, more 
high-resolution studies of K-Pg boundary records around the globe 
are needed to confirm the degree of hemispheric heterogeneity in 
the post-impact biotic recovery. To conclude, our results highlight 
that the photosynthetic shut-down induced by the large volume of 
silicate dust with grain sizes between ~0.8 and 8.0 um, together with 
additional effects of sulfur and soot, probably led to a disastrous col- 
lapse of primary productivity in land and ocean realms, steering the 
global mass extinction at the K-Pg boundary. 


Online content 

Any methods, additional references, Nature Portfolio reporting sum- 
maries, source data, extended data, supplementary information, 
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Methods 

Laser-diffraction grain-size analysis of K-Pg boundary 
sediments 

Laser-diffraction grain-size analysis was performed on 40 sediment 
samples from the X-2741-A and X-2761 (Supplementary Table 1) sec- 
tions from the Tanis K-Pg site (North Dakota, United States”) to deci- 
pher the depositional processes of the Upper Cretaceous point bar 
strata, unit K, the event deposit, units 1 and 2, as well as the Palaeocene 
intervals, P1 and P3 (Extended Data Fig. 1b). The analysis took place at 
the Sediment Laboratory at the Vrije Universiteit Amsterdam using a 
Sympatec HELOS KR laser-diffraction particle-size analyser’. Before 
particle-size analysis, the sediment samples were prepared following 
a procedure” to ensure solely the lithogenic fraction was measured. 
Organic matter and carbonate material were removed by treating 
the samples with 5-10 ml 30% H,O, and 5-10 ml 10% HCI, respec- 
tively, both followed by heating to boiling point. After decarboniza- 
tion, the Palaeocene lignite samples (unit P2) did not yield sufficient 
lithogenic fraction to achieve accurate results. To prevent clumping 
of clay minerals to agglomerates, 300 mg sodium pyrophosphate 
(Na,P,0,-10H,O) was added to the suspension, followed again by 
heating until boiling point. After cooling, the samples in suspension 
were measured in the HELOS laser-diffraction sensor device, with 
the resulting lithogenic grain-size distribution data (ranging from 
0.1 to 2,000 um, measured in 56 classes) expressed in volumetric 
percentage. 

Aclear two-step fining upward sequence is recognized in the Tanis 
K-Pg boundary succession (Extended Data Fig. 1b), which has been 
attributed to the arrival of two successive tsunami-like surge waves 
generated by the Chicxulub impact event“. The sample just below the 
Palaeocene lignite, together with the other two K-Pg claystone samples, 
reveals a distinct and uniform particle-size distribution with ~80% 
clay fraction, indicative of fine atmospheric settling”. This sharply 
contrasts with the Tanis event deposit, which shows large heterogenei- 
ties with more pronounced bimodal distributions, including a coarser 
peak of fine sand around 125 um, indicative of transport by water” 
(Extended Data Fig. 1c). The uppermost K-Pg claystone sample yields 
avery fine median of 2.88 um (Fig. 2a and Extended Data Fig. 1c), and 
this grain-size distribution was used in our numerical modelling study. 
Note that the measured mass density spectrum is fitted by a trimodal 
log-normal size distribution (see the model parameters of each mode; 
Extended Data Fig. 2a). 


Mass to number density spectrum conversion 

This analysis transforms the mass density spectrum of K-Pg grain-size 
data into the number density spectrum (Extended Data Fig. 2), which 
provides a compatible input for a GCM”. Let Q,(x;) be a cumulative 
particle-size distribution from a given sediment sample: 


AQ, Xi-1: Xi) = Q; Xi) — Qr Xi-1) (1) 


In equation (1), x; is the particle size and AQ,(x,) is a differential 
particle-size distribution. The spectrum density, g,(x)), is given by: 


AQ, (x) _ Q, (x) = Q, Xi) 


AX; Xj —Xi-1 


(2) 


Or (Xj) = 


Herer denotes the type of density spectrum: r= 0 (number) or r=3 
(mass). To convert q,(x;) fromr= 3 tor = 0, we followed a well-described 
method”. First, a given density spectrum q,(x,) can be transformed 
into g,(x;) by: 


qr x) = X qe (xj) IM, —e,e (3) 


In equation (3), M,—e isthe moment of distribution. The moment 
of particle-size distribution can be generalized by: 


Xi,max 


Meesr— 
Mer = xkqr x) dx; = ee (4) 


r—e,e 
Xi,min 


Here My, is the Ath moment of r-type density spectrum, q,(x;). 
Since the transform will be applied from mass to number spectrum, 
recalling equation (3) and taking e= 3 and r=0: 


qo (Xi) = X7 q3 Xi) /M_33 (5) 


Supposing alog-normal particle-size distribution, the kth moment 
of a spectrum density, g,(x;), can be expressed as follows: 


k an) 
Mkr = Xgr EXP (Sm osr) (6) 


In equation (6), w,and o, denote the median and standard devia- 
tion of particle-size distribution, respectively. Substituting k = -3 and 
r=3into equation (6), we have the following expression: 


2 


M_33 = X;3 exp ( C3) 


2 E 2 
za 3 In oe = x33 exp (4.5In oga) (7) 


Then substituting equation (7) into equation (5), the final form of 
transformation is as follows: 


qo (Xi) = x7? B) 


(8) 


as exp (4.5In°o,;) 


In equation (8), q(x) isthe mass density spectrum and qo(x;) isthe 
number density spectrum. 


Modelling the log-normal number density spectrum 

This was carried out by using the output of laser-diffraction analysis, 
gx), referring to the mass density spectrum on the logarithmic 
abscissa (Extended Data Fig. 2a). Here, 3 (x,)is defined by the logarithm 
of base 10 that can be expressed as: 


fig (log x;) = 43 i) (9) 


By converting equation (9) into the natural logarithm: 


AG (Inx;) dInx; = ā% (log x;) dlogx; > AG (Inx;) = ñ (log. x;) /2.303 (10) 


Letting D, =x,is more convenient to denote the particle size. Then 
fig(InD,), the mass density spectrum, can be transformed to number 


density spectrum, Ā%(InD,); by substituting into equation (8), we have 
following: 


ñg (IND 
ii, (InD,) = D5? — My (InDp) 
D7 g3 EXP (4.5ln? 0,3) 


(11) 


To satisfy a probability density function, equation (11) can be 
normalized by: 


fy (InD,) 


ny (InDy) = —s (12) 
Ay, (InD,) dinD, 


=% 


fiy(InD, ) The term can now be approximated by alog-normal probability 
density function. 
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In equation (13), Nrefers tothe number concentration of particle- 
size distribution. The exact expression needs to be given in the form 
of n,(D,), instead of ñ% (InD,). 


dD, 
ny (InD,) dinD, = ny (Dp) dD, > nj, (InD,) = ny (Dp) dnb, = Ny (Dp) Dp 
(14) 


The term ny (D,)in equation (14) reads as follows: 


nê (InD,) N 
Ny (Dp) = N D P exp 
Pp V2nD,Inog,o 


2 
(InD, — InDpg.o) (15) 
2In? dg,9 


Equation (15) can be rewritten in differential form. 


2 
(InD, — InDpgo) a6) 
2In’ og. 


Ny (Dp) dD, 


: exp ( 
V2nD,|ndg, 


To satisfy alog-normal probability density function, equation (16) 
is scaled by a weighting factor, wo. 


2 
ny (Dp) dD, = (rip lagn) ) (17) 


2 
2In" ogo 


Wo ( 
exp 
V2nD, Inozo 


Concisely, the procedure we followed is that the mass density 
spectrum from the laser-diffraction grain-size analysis, q3(x;), is first 
converted to a number density spectrum via equation (11), which is 
in turn modelled through alog-normal probability density function 
using equation (17), to use as an input in GCM simulations. The model 
parameters in equation (17) were determined as follows (Extended 
Data Fig. 2b): 


Wo = 7.24, Dpg,o = 0.125 um, Inozo = 0.446 (18) 


Here, D,zo and Ing, are the median and logarithmic standard 
deviation of modelled number density spectrum, respectively. 


GCM overview 

We performed palaeoclimate GCM simulations by our paleoEarth 
implementation of the general-purpose planetWRF model”. The plan- 
etWRF model has been widely applied and validated for the modelling 
of planetary atmospheres**” on the basis of the extensively used 
terrestrial atmospheric model, Weather, Research and Forecasting 
(WRF) model”. Here our paleoEarth model” is a GCM that is special- 
ized to simulate the Chicxulub impact at the K-Pg boundary, using the 
latest Cretaceous palaeography for geographic boundary conditions”. 
Our GCM makes use of a simplified plant functional type scheme for 
land-surface conditions, approximating the latest Cretaceous veg- 
etation”. It resolves the atmospheric transport of fine-grained ejecta 
(consisting of sulfur, soot and silicate dust) generated by the Chicxu- 
lub impactor, considering the radiative and microphysical feedback 
of climate-active ejections. Regarding the verification of GCM, the 
simulated pre-impact surface temperatures agree with proxy-based 
latest Cretaceous temperature reconstructions (Extended Data Fig. 4). 
Moreover, inarecent intercomparison study”, our initial GCM recon- 
structions for the post-impact global temperature drop fall in line with 
those of previous GCM simulations following the Chicxulub impact. It 
is worth mentioning that our initial GCM” is further improved here. In 
the initial GCM”, a simple algorithm is used for the dry deposition of 
tracers; however, in the present version, we used an advanced atmos- 
pheric advection scheme of tracers available in the planetWRF. We 
updated the optical properties of soot and dust, such as the specific 


extinction, single-scattering albedo and asymmetry factor as a func- 
tion of wavelengths, which led to less surface cooling with regard to 
the climatic response of dust and soot emissions. Moreover, soot par- 
ticles were injected only from the impact site, up to the stratosphere 
top in our initial GCM”, which we updated here in the present version 
injecting soot only over the land surfaces from the surface up to lower 
stratosphere (‘Emission configurations’). Note that we added an algo- 
rithm into the present GCM that conserves the global mass of emitted 
tracers (sulfur, soot, or dust) over the whole GCM domain from one 
time step to the next one. Consequently, the global mass of tracers is 
precisely conserved in the atmosphere, sustaining monotonic tracer 
deposition. The paleoEarth GCM consists of three main components 
besides the dynamical core of planetWRF. The first component is a 
module implemented for the asteroid/cometary impact dynamics, 
in which the transport of impact-generated fine-grained ejecta (con- 
sisting of sulfur, soot and silicate dust) is governed in the GCM. The 
second and third components are radiative transfer and microphysics 
modules, respectively. 


Radiative transfer and microphysics modelling 
The radiative transfer module, the second component, is used to 
include the radiative effect of impact ejecta, such as the absorption, 
reflection and scattering by ejected particles in the short-wave and 
long-wave spectra. It was built within the NASA (National Aeronautics 
and Space Administration) Goddard short-wave/long-wave radia- 
tion scheme®*, that is, composed of eight spectral bands within the 
ultraviolet and visible range (<700 nm) and three spectral bands in 
the infrared region (>700 nm). In this component, specific extinc- 
tion, single-scattering albedo and asymmetry factor of sulfate were 
set as functions of wavelengths in short-wave and long-wave spectral 
bands®. Regarding the dust, optical properties are specified following 
previous studies in the short-wave” and long-wave” spectrum. Optical 
properties of soot are determined in short-wave and long-wave spectral 
bands*. The microphysics module, the third component, resolves the 
microphysical processes of fine-grained ejecta, such as dry and wet 
deposition processes of soot, sulfur and silicate dust. Microphysical 
processes are crucial in determining deposition rates and, as a result, 
the lifetimes of deposited particles in the atmosphere. For the dry 
deposition modelling of soot and dust, our GCM considers relative 
contributions of gravitational and aerodynamic forcing as well as the 
Brownian diffusion, impaction and interception of particles, using a 
particle-size-aware resistance model®”. 
Va = Vg + (Ra + Rs + Ra Rs Vg)” (19) 

Here, V,and V, are the deposition and gravitational settling veloci- 
ties, while R, and R, are aerodynamic and surface deposition resist- 
ances. The latter term reads as follows: 


1 
€oU,.(Ep + Eim + Ein) 


(20) 


Rs 


In equation (20), £= 3 is an empirical constant®’ and u. is the 
friction velocity of near-surface winds that the GCM computes. Here 
E,, Eim and En are collection efficiencies related to Brownian diffu- 
sion, impaction and interception of particles determined from previ- 
ous models”. The dry deposition rate of sulfate on land was set to 
0.1 cm s”, referring to the value of coniferous and deciduous forest 
covering land”. Likewise, the ocean deposition rate is prescribed to be 
0.1 cm s”, which falls between the 0.05 cm s~” (ref. 73) and 0.20 cms? 
(ref. 65) reported in previous studies. Note that our sulfur-cycle simu- 
lations do not include photochemistry reactions. We presumed that 
the sulfur injection is fully converted into sulfate as the maximum 
scenario (‘Emission configurations’) as sulfur has around three to four 
times higher deposition rates on land and inthe ocean, 0.6-0.8 cms? 
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(ref. 65), leading to faster deposition than sulfate. Moreover, the wet 
deposition of particles is modelled by resolving in-cloud rainout and 
washout processes” ”°. For the dust microphysical modelling, our GCM 
takes dust lifting into account, consisting mainly of three mechanisms: 
(1) dust lifting by aerodynamic forces, (2) saltation bombardment 
and (3) breakdown of large aggregates (disaggregation)’””’. Note that 
the microphysical modelling of soot and dust were implemented in 
a size-resolved way” unlike the monodisperse assumption in sulfate 
modelling. Here the size-resolved modelling follows both our recent 
two-moment dust transport scheme initially derived for the Martian 
dust cycle“ and atwo-moment framework widely used in community 
atmosphere models”. Our soot and dust microphysics models do not 
include coagulation microphysics, which might pose as important 
for the transport of aerosols in the atmosphere. Nanometric particles 
below 0.1 um, that is, 0.015 um < D, < 0.052 um, whose range is referred 
to as the Aitken mode”, are formed by two processes: (1) condensa- 
tional growth on existing aerosol particles and (2) coagulation due 
to the random particle collisions. These nanometric particles can 
further grow into larger particles or chains, resulting in the so-called 
accumulation mode (0.056 um < D, < 0.26 pm)”, in which coagulation 
can still occur especially at high particle concentrations following 
the K-Pg impact. The median diameters of dust (0.125 um) and soot 
(0.22 um) in our simulations are larger than the Aitken-mode interval 
while lying within the range of accumulation mode (Extended Data 
Fig. 3). In our GCM, the inclusion of a coagulation scheme would have 
formed larger soot aggregates up to 0.26 um (the upper range of accu- 
mulation mode). Ref. 23 put forward that the coagulation can increase 
short-wave optical properties owing to the large soot emissions. Such 
particle sizes (up to 0.26 um) are close to the median diameter of soot 
(0.22 um), falling well into the high-probability range of log-normal 
size distribution used in our size-aware two-moment aerosol scheme. 
Therefore, any alteration in particle deposition and atmospheric life- 
time or radiative response would be minor. This may also apply to the 
dust case. Even if larger dust aggregates with the same upper range (up 
to 0.26 pm) would form by the coagulation, the median size of dust 
particles (0.125 um) would be close to the given size range. Despite a 
slight underestimation in current radiative response, such a change in 
particle size up to 0.26 pm would not have substantial influence onthe 
atmospheric settling processes on land and ocean (Fig. 4) and thus on 
our main findings. Therefore, we assume the effect of coagulation to 
be lesser on the overall atmospheric processes of the given soot and 
silicate dust ejecta following the Chicxulub impact. 


Description of simulation set-up and model physics 

For the land-surface physics, in our GCM simulations we made use of 
the five-layer thermal diffusion scheme*®’. To model the water-cycle, 
cloudand precipitation microphysics, we used the Purdue-Lin micro- 
physics scheme“. The modified Tiedtke scheme is utilized for the 
cumulus parameterization®***’. The planetary boundary-layer turbu- 
lence is modelled by means of our recent planetary boundary-layer 
scheme™. The atmospheric surface layer is modelled by using the 
revised MMS scheme®. Regarding ocean modelling, we use the Pollard 
ocean mixed layer model**” that accounts for the wind-driven turbu- 
lent mixing and deepening in the ocean. As we are interested mainly 
inthe relative roles of different aerosols, we do not expect that using a 
more-sophisticated ocean and/or land surface would change our main 
conclusions. Moreover, we imposed the latest Cretaceous climatic 
conditions, taking global-average atmospheric CO, concentration to 
be 560 parts per million™'®. In our GCM simulations, we did not release 
CO, due to impact heating and global wildfires. It has no dominant 
effects on short timescales as focused in the present study, leading to 
along-term warming of about 100,000 years since the impact winter". 
Although it could induce global warming shortly after the K-Pg impact, 
the extent of warming is minor. Simulations” indicate that the change 
of C release from 115 to 1,615 Gt can result in -2 °C warming from 5 to 


20 years after impact. For the orbital forcing, a circular orbit with an 
obliquity of 23.5° and a solar constant of -1,354 W m”? are assumed”, 
similar to other palaeoclimate studies. The moment of the Chicxu- 
lub impact, thus the release of fine-grained ejecta, is assumed to be 
initialized in the boreal spring”. This assumption is based on recent 
osteohistological and isotopic studies from uniquely preserved fossil 
fish from the Tanis K-Pg site, and this moment of impact probably sub- 
stantially influenced selective biotic survival across the K-Pg boundary 
catastrophe’. The horizontal model resolution of GCM is 5° x 5° over 
longitudinal and latitudinal directions, having 27 vertical sigma layers 
extending through the stratopause. Furthermore, the time integration 
for eachimpact simulation is carried out for 35 years following an initial 
spin-up simulation of 15 years, in which the latest Cretaceous condi- 
tions stabilized (Extended Data Fig. 10). 


Calculation of PAR flux 

PAR flux refers to the solar radiation reaching the surface intheA = 400- 
700 nm spectral range, indicating how intense the radiative flux is to 
sustain photosynthesis being used by photosynthetic organisms**. It is 
thus highly important in quantifying the global primary production™. 
In our paleoEarth GCM, the time evolution of global PAR flux is com- 
puted within the radiative transfer module at each model time step. As 
the NASA Goddard short-wave/long-wave radiation scheme®* solves 
the radiative transfer in spectral bands, the eighth spectral band corre- 
sponds tothe PAR spectral range, in which we saved model results at each 
hour and post-processed for detailed interpretations (Figs. 3d and 4). 
Our simulations show that the global PAR flux is 89 W m7 in the latest 
Cretaceous, which is larger than the present-day fluxes (66-74 W m°). 
Here we computed the present-day range following the approach” in 
which the PAR flux is converted as the fraction (k, = 0.45-0.5) of global 
solar irradiance reaching Earth’s surface (R,). The annual mean of 
global-average R, was taken to be 147.3 W m”? from the ten-year radia- 
tion budget of Earth between 1984 and 1993 (ref. 91). 


Emission configurations 

The particles were emitted instantaneously at the moment of impact. 
In our GCM, the latest Cretaceous conditions stabilize following an 
initial spin-up simulation of 15 years. The next year (t = 16 yr) hosts the 
impact event when the day of year is 151, falling into the boreal spring 
season’*®. Ref. 10 reported that the fast-moving ejecta cloud could carry 
dust, soot and sulfate particles across Earth within four to five hours of 
impact. Silicate dust and sulfur spread are slower due to atmospheric 
transport; nevertheless, our GCM results show almost a global cover- 
age already at one day after impact (Fig. 5). Here we injected sulfur 
and silicate dust from the impact site uniformly in the atmospheric 
column and up tothe model top, that is, top of the stratosphere, while 
soot particles were injected over the whole land surfaces following 
ref. 18. The current emissions design of silicate dust and sulfur may 
affect aerosol concentrations within the first few days after impact 
(up to the first week), overestimating concentrations near the impact 
site compared with distant locations. Enormous quantities of sulfur 
(~325 Gt) and silicate dust (~2,000 Gt) travelling at hypervelocities 
exceeding 5 km s” (ref. 10) could suggest that the silicate dust and sulfur 
can rather be ejected globally asin ref. 39. Nevertheless, emitting from 
the impact site instead of injecting globally would not have a notable 
influence on the timescale focused in our study (from first few weeks 
to 25 years). Concerning the emissions design of soot, for the vertical 
distribution, ref. 39 suggested a Gaussian vertical profile consisting of 
two peaks centred at the tropopause and surface. This way, the normal- 
ized injection rate at the surface (lower peak) is nearly three times larger 
than that at the tropopause (upper peak), in which the upper peak has 
a half width of 3 km injecting up to the altitude of 25 km. We injected 
soot particles from the surface up to 25 km (lower stratosphere), while 
assuming a uniform profile. Our approach is based onthe assumption 
that the rapidly moving ejecta cloud, comprising large quantities of 
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soot (~20.8 Gt ref. 18)) from fires, disperse quickly in the vertical lay- 
ers of atmosphere soon after impact. Strong vertical motions at high 
altitudes as well as up- and downdraughts, entrainment and turbulence 
inthe boundary layer” could carry the smoke likely well mixed in tropo- 
sphere and lower stratosphere. More-sophisticated injection profiles of 
sootas in ref. 39 could better estimate the atmospheric perturbations 
locally and vertically in the immediate aftermath of impact (first days 
to the first two weeks). However, it may have no substantial effect on 
the longer timescales, which we focused on in this study. Regarding the 
sulfur emission, given the relatively high rates of dry and wet deposi- 
tion, therefore the reduced atmospheric lifetime of SO,, we assumed 
that SO, is fully oxidized to form sulfate aerosols as the upper limit of 
atmospheric perturbations. The atmospheric lifetime of sulfate can 
be longer than our upper estimate as the conversion rate from SO, 
to sulfate may take a while" due to homogeneous or heterogeneous 
chemical pathways”. However, the water release from the shallow sea 
impact site” can speed up this conversion rate”. Moreover, a fraction 
of SO, can convert into more quickly removed SO, inthe atmosphere, 
increasing the SO./SO, ratios within the impact vapour cloud. Impact 
experiments demonstrated that the major portion of sulfur in the 
impact ejecta occurs in the form of SO, (ref. 96). If this is the case, it 
implies that the conversion to sulfate aerosols would form even more 
rapidly in the stratosphere”. Note that the vertical extent of our GCM 
spans fromthe planetary surface to the top of the stratopause. Here we 
focused onthe relative roles of different impact emissions on the over- 
all atmosphere with coupled troposphere-stratosphere. Individual 
responses of troposphere and stratosphere were not investigated in 
the present study; however, the stratospheric aerosols would probably 
have a longer atmospheric lifetime than the tropospheric particles, 
which settle out quickly”. Since some fraction of impact-induced 
aerosols can reach higher altitudes and be influenced through photo- 
chemical processes, the stratosphere and mesosphere/thermosphere 
coupling following the aftermath of the Chicxulub impact event shall 
be addressed by further studies. 


Data availability 

The palaeoclimate general circulation model (GCM) output data as well 
as the silicate dust grain-size data from the Tanis K-Pg site that support 
the findings of this study are publicly available in the OSF repository 
via https://doi.org/10.17605/OSF.IO/2CDQG. The proxy-based latest 
Cretaceous temperature reconstruction data are available online in 
the PANGEA repository: https://doi.org/10.1594/PANGAEA.879763. 


Code availability 

The Python and Matlab source codes developed for reproducing the 
figures in this study are publicly available at the GitHub repository via 
github.com/cem-berk-senel/naturegeoscience-chicxulub/. The Plan- 
etWRF model is available upon request from https://planetwrf.com/. 
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paleogeographic reconstruction and relief map for the latest Cretaceous’ as 
used in this modeling study with locations of Chicxulub and Tanis indicated. 
The base map is based on the latest Cretaceous paleogeographic data”. 
b Stratigraphy of the Tanis K-Pg boundary event deposit highlighting the 
lithological units (adapted from“; based on sections X-2741-A and X-2761) 
together with data on grain-size classes (clay, silt, and sand fractions), median 
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um (Toon et al., 2016”; Wolbach et al., 1985”). Concerning the aerosol life cycle 
and processes, coagulation is one of the crucial microphysical mechanisms 

that might pose as important for the transport of aerosols in the atmosphere. 
Nanometric particles below 0.1 um, that is, 0.015 um < D, < 0.052 um”, whose 
range is referred to as the Aitken-mode (within the cyan dashed lines), are formed 
by two processes: (i) condensational growth on existing aerosol particles, 


and (ii) coagulation due to the random particle collisions. These nanometric 
particles can further grow into larger particles or chains, resulting in the so-called 
accumulation-mode (0.056 um < D, < 0.26 um)” (within the red dashed lines) 
where the coagulation can occur especially at high particle concentrations 
following the K-Pg impact. The median diameter of soot (0.22 um”°””) in our 
simulations are prominently larger than the Aitken-mode interval, while lying 
within the range of accumulation-mode. 
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lines display the zonal mean of land temperatures during the boreal summer and 
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a Latest Cretaceous: 1y before impact (annual mean) d Impact winter: 6m after impact 


b Latest Cretaceous: 1w before impact (boreal spring) e Impact winter: 2y after impact 


c Impact winter: 1m after impact f Impact winter: 10y after impact (annual-mean) 
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Extended Data Fig. 5 | Global surface temperature reconstructions using b Latest Cretaceous, 1 week before impact. c Impact winter, 1 month after impact. 
the combined fine-grained ejecta scenario. Results are displayed ona latest d 6 months after impact. e 2 years after impact. f10 years after impact (annual- 


Cretaceous paleogeographic map (Extended Data Fig. 1a) and shown for different mean). Base maps are based on the latest Cretaceous paleogeographic data”. 
time snapshots. a Latest Cretaceous, 1 year before impact (annual mean). 
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Extended Data Fig. 6 | Impact-generated global surface net radiative The solid purple dashed line denotes the moment of Chicxulub impact, thatis, 
responses. The temporal evolution from -2 years of the latest Cretaceous boreal spring season”. Our paleoclimate simulations indicate that the drastic 
towards 25 years of post-impact conditions, for the individual silicate dust, changes in surface net shortwave/longwave radiation stabilize to pre-impact 
sulfur, soot, and combined scenarios. a Global-average surface net shortwave levels within the first 3 years after impact. Accordingly, this timescale, in which 
radiation flux. b Global-average surface net longwave radiation flux. Here in large radiative anomalies emerged, determines the timescale of the initial 
x-axis, the year of O refers to the start of the year where the impact event occurs. extreme cold (Fig. 4a). 
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a Latest Cretaceous: 1 day before impact (boreal spring) 


Global photosynthetic active radiation [W/m?] 
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purple colorbar represents the photosynthetically high and low radiative flux, 
varying between 0-160 W/m”. Base maps are based on the latest Cretaceous 


Extended Data Fig. 7 | Global PAR flux reconstructions in the latest 
Cretaceous. Land-Ocean PAR flux a1 day before impact (boreal spring) and 
b6 months before impact (austral summer), displayed ona latest Cretaceous paleogeographic data”. 
paleogeographic map (Extended Data Fig. 1a). The range of the green and 
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Global photosynthetic active radiation [W/m?] 
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Extended Data Fig. 8 | PAR flux reconstructions following the Chicxulub onalatest Cretaceous paleogeographic map (Extended Data Fig. la). The range of 
impact. Land and ocean PAR flux from 1 day (post-impact state, instantaneous)to1 the green-white colorbar denotes the photosynthetically high and low radiative 
and 2 weeks after impact for a silicate dust; b sulfur; and c soot scenarios, displayed flux. Base maps are based on the latest Cretaceous paleogeographic data*’. 


Nature Geoscience 


Article 


https://doi.org/10.1038/s41561-023-01290-4 


Global column-integrated fine-grained ejecta mass 


1016 T T 


1014 
1012 
1010 


Chicxulub impact 
u> 


[kg] 
= 
D 


10? pae 
10° Latest 


2 Cretaceous 


oe 


= Tanis K-Pg silicate dust (this study) 

= Dust nanoparticles (Toon et al. 2016) 

— Dust type 2 coarse spherules (Toon et al. 2016) 
= = Shocked K-Pg quartz grains (Pope et al. 2002) 
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Extended Data Fig. 9 | Effect of silicate dust particle size on the global 
column-integrated fine-grained ejecta mass. Orange line refers to the present 
study, using the Tanis K-Pg silicate dust. Gray and black lines show GCM results 
using particle size constraints reported in previous modeling studies'*”’. The 
gray line displays the response of nanometric sized particles, indicating a 
deposition rate with an atmospheric lifetime of -7 years (as a lower threshold). 
The black line refers to the type 2 spherules”, representing microkrystites of 
250 um in diameter prone to very swift gravitational settling within a few days 
after impact. Cyan dashed line displays the shocked ejected quartz grains (mean 
diameter of 50 um) defined as clastic debris”. We use the same amount of ejecta 
release in each GCM simulation, in the order of 2x10" g as an upper limit. The 
optical properties of nanoparticles and type 2 spherules are the same as inthe 
Tanis K-Pg silicate dust simulation, as we compare the microphysical response to 
the changes in particle size. Here in x-axis, the year of O refers to the start of the 
year where the impact event occurs. The purple dashed line denotes the moment 


of Chicxulub impact, that is, boreal spring season”. Regarding nanoparticles, 
those nanometric sized particles (median diameter of 20 nm) would growinto 
larger particles in atmosphere due to the coagulation. Such larger aggregates 
would have lower deposition rates on land and ocean (Fig. 3), hence higher 
atmospheric lifetimes. To illustrate, the deposition rate of nanoparticles (gray 
dashed line) would have occasionally shifted rightward through the response 

of silicate dust (orange dashed line) depending on the rate of coagulation. 
Therefore, the present simulation of nanoparticles, excluding coagulation, 
would serve as the minimum threshold for the atmospheric lifetime (t - 7 years). 
The inclusion of coagulation mechanism forming larger aggregates would 

lead to lower deposition rates on land and ocean (Fig. 3) for some fraction of 
nanoparticles, thus relatively high atmospheric lifetime of more than 7 years. 
Nevertheless, we do not expect nanoparticles to have an atmospheric lifetime 
and PAR response as substantial as single soot or micrometer-sized silicate dust. 
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Extended Data Fig. 10 | Global-average surface temperature. It is same as Fig. 4a, yet the time evolution is shown from 15 years before the Chicxulub impact instead 
of 2 years, for the individual silicate dust, sulfur, soot, and combined scenarios. The first 15 years correspond to the model initial spin-up simulation of 15-years, 
in which the latest Cretaceous conditions stabilized. 
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